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ABSTRACT 

The spectral energy distributions (SEDs) of dust-enshrouded galaxies with powerful restframe far- 
infrared(IR) emission have been constrained by a range of ground-based and space-borne surveys. 
The IRAS catalog provides a reasonably complete picture of the dust emission from nearby galaxies 
(redshift z ~ 0.1) that are typically less luminous than about 10 12 L@. However, at higher redshifts, 
the observational coverage from all existing far-IR and submillimcter surveys is much less complete. 
Here we investigate the SEDs of a new sample of high-rcdshift submillimcter-selected galaxies (SMGs) , 
for which redshifts are known, allowing us to estimate reliable luminosities and characteristic dust 
temperatures. We demonstrate that a wide range of SEDs is present in the population, and that a 
substantial number of luminous dusty galaxies with hotter dust temperatures could exist at similar 
redshifts {z ~ 2 — 3), but remain undetected in existing submillimeter surveys. These hotter galaxies 
could be responsible for about a third of the extragalactic IR background radiation at a wavelength of 
about 100 /im. The brightest of these galaxies would have far-IR luminosities of the order of 10 13 L© 
and dust temperatures of the order of 60 K. Galaxies up to an order of magnitude less luminous with 
similar SEDs will be easy to detect and identify in the deepest Spitzer Space Telescope observations 
of extragalactic fields at 24 /im. 

Subject headings: galaxies: evolution - methods: observational - galaxies: interstellar medium - 
galaxies: starburst - cosmology: observations 



1. INTRODUCTION 

Deep submillimcter-wavc surveys for distant galaxies 
offer direct access to high redshifts, because their sensi- 
tivity is almost independent of galaxy luminosity for a 
fixed spectral energy distribution (SED) over the wide 
redshift range 1 < z < 5 (Smail et al. 2002). The 
positional accuracy of these surveys is only ~ 8 arcsec, 
and optical counterparts are typically faint. Hence it 
has been difficult to identify and measure the redshifts 
of a substantial fraction of these submillimeter-selected 
galaxies (SMGs). 

These difficulties have now been overcome with rest- 
UV spectroscopy of SMGs located in very deep (~ 
10/iJy) radio images by Chapman et al. (2003b; 2004a). 
Redshifts have been found for about half of the SMC 
population using the Keck I telescope's blue-sensitive 
multiobject spectrograph LRIS-B. The resulting sample 
of 73 spectroscopic redshifts is very valuable, allowing a 
fairly accurate restframe SED to be determined for each 
galaxy, under reasonable assumptions about their radio- 
far-IR properties. The SED model is most strongly af- 
fected by the value of the effective dust temperature of 
the galaxies Td, which can be treated as a proxy for the 
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peak restframe wavelength of the SED. Just two parame- 
ters, Td and the far-IR luminosity L? allow the parameter 
space of SEDs that describe observational data at wave- 
lengths longer than about 20 /im to be spanned reason- 
ably well (Dale et al. 2001; Blain, Barnard & Chapman 
2003; hereafter BBC03). Two other parameters are re- 
quired to describe the detailed shape of the SED. The 
Rayleigh- Jeans spectral slope of dust emission is fixed 
by the emissivity index of dust f3 (/„ oc v 2+l3 ). The 
value of (3 is degenerate with temperature, in the sense 
that models in which Td/3 is approximately constant pro- 
vide similarly accurate descriptions of the SED (see Fig. 2 
of BBC03). This effect could be important if Td were 
treated as a real physical temperature, but when Td is 
used as the key descriptive parameter in the SED the 
effect of this degeneracy between the fitted values of Td 
and P is not important. The slope of the mid-IR spec- 
trum /„ cx v a is set by a second parameter, which has 
little effect on the value of L at fixed Td. 

Spectroscopic redshifts are especially important for un- 
derstanding the SEDs of SMGs as there is a very signif- 
icant, pernicious degeneracy between the inferred dust 
temperature Td and redshift z, which only allows ac- 
curate measurement of the quantity Td/(1 + z) (Blain 

2 Defined as the luminosity between wavelengths of 3 mm and 
10 /im, and typically close to the bolometric luminosity. 
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1999; BBC03). The large Chapman et al. sample of se- 
cure spectroscopic redshifts for SMGs provides a way to 
break this degeneracy, and thus unambiguously reveals 
the range of SED properties of the high-redshift popula- 
tion of SMGs. This is the population of galaxies for which 
radio-submillimeter-far-IR photometric redshifts are of- 
ten sought (Carilli & Yun 1998; Hughes et al. 2002; Yun 
& Carilli 2002; Aretxaga et al. 2003; Wiklind 2003), and 
we use this new, accurate SED information, to address 
the precision of these radio-far-IR photometric redshift 
techniques. We then compare the SEDs of low- and high- 
redshift samples of luminous dusty galaxies to highlight 
a region of the L-T^ parameter space where unrecog- 
nized populations could still lurk, and discuss ways to 
find them. 

Throughout this paper we assume a cosmology with 
O = 0.3, tt A = 0.7 and H a = lOO/ikms" 1 Mpc" 1 with 
h = 0.65 to remain consistent with earlier work. In the 
WMAP cosmology with Q, = 0.27, fl A = 0.73 and h = 
0.71 the quoted luminosities would be reduced by a factor 
of about 15% at z ~ 2.5. 

2. LUMINOSITIES AND TEMPERATURES OF 
DUST-ENSHROUDED GALAXIES 

In Fig. 1 we show values of inferred dust tempera- 
ture Td and bolometric far-IR luminosity L for sev- 
eral samples of dust-enshrouded galaxies. These include 
well-studied low-redshift IRAS galaxies with submillimc- 
ter measurements (Dunne et al. 2000; Dunne & Eales 
2001), moderate-redshift galaxies selected jointly from 
the IRAS and the VLA-FIRST survey (Stanford et al. 
2000) and high-redshift galaxies selected using both In- 
frared Space Observatory (ISO) at 15 fim and radio im- 
ages (Garrett et al. 2001). Several examples of serendip- 
itously selected high-redshift SMGs with known redshifts 
are also shown, including some gravitationally lensed ex- 
amples. The values of T<j and L inferred from these sam- 
ples and their uncertainties were discussed by BBC03, 
and have now been updated for Fig. 1. We now also plot 
the inferred Td and L values derived for the Chapman et 
al. (2003b, 2004a) sample of 73 SMGs with spectroscopic 
redshifts. 

Following BBC03, values of Td and L are derived as- 
suming a modified blackbody SED with a single dust 
temperature, and a dust emissivity index (3 = 1.5, which 
is extrapolated with a continuous gradient to a power- 
law f v oc v a at frequencies above which the gradient of 
the blackbody spectrum is steeper than this power law. 
The value of a. = —1.7 is used here, because it provides 
a good description of the shape of 15-/um source counts 
based on the form of evolution inferred for the SMGs 
(Blain et al. 2002). More information about the range 
of values of this parameter awaits observations with the 
Spitzer Space Telescope. 3 Arp 220 has amongst the steep- 
est indices known with a = —2.9 (BBC03). 

All the galaxy samples shown in Fig. 1 have been fitted 
using the same function, so the derived luminosities and 
temperatures for each sample should be directly com- 
parable. This three-parameter (T^i, a, (3) model is the 
simplest that can describe the data adequately, and Td is 
the most important parameter (BBC03). To emphasize, 
the value of Td does not necessarily reflect the true phys- 



ical temperature of most of the dust in the galaxy, but 
does provide a good description of the general properties 
of the SED. 

When fitting the SEDs, we assume that the low- 
redshift radio-far-IR flux density correlation (Condon 
1992), resulting from the fact that high-mass stars are 
the common power source for both synchrotron radio 
emission via supernovae shock acceleration, and far-IR 
dust emission via absorbed blue/UV stellar light, applies 
to the Chapman et al. (2003b, 2004a) galaxies. Un- 
der this assumption, we can determine a temperature 
Td to describe the SED from just two data points at 
1.4 GHz and 850 /im. If there are systematic changes 
in the radio-far-IR flux density ratio or its scatter with 
redshift, then these values will be in error. So far there 
is no evidence for evolution in the form of the far-IR- 
radio correlation (Garrett et al. 2001). Information to 
test the inferred SEDs could be provided by difficult-to- 
obtain shorter submillimctcr flux density measurements 
in the 350- and 450-^m atmospheric windows, or from 
space-borne far-IR data. So far, archival 450-yum SCUBA 
measurements at the James CLerk Maxwell Telescope 
(Smail et al. 2002) and new 350-/xm measurements using 
the Caltech Submillimctcr Observatory (Mauna Kea)'s 
SHARC-2 camera by Kovacs et al. (2004) for about 
10 SMGs show no major discrepancy as compared with 
the radio-submillimeter result: the dust temperatures in- 
ferred from the measured submillimeter color between 
350 and 850 /im are consistent with those derived using 
the far-IR-radio correlation (Fig. 1) to within the errors. 

The presence of an active galactic nucleus (AGN) 
might be expected to boost the radio flux density, leading 
to an overestimate of Td (Yun, Reddy & Condon 2001). 
This is clearly the case in only 3 out of 73 galaxies in 
this sample, for which very high dust temperatures and 
luminosities are inferred. Where X-ray observations are 
available for the Chapman et al. sample, the results are 
inconsistent with a large fraction of Compton-thin AGNs 
(Alexander et al. 2004). 

Note that, as discussed by BBC03, the addition of more 
continuum data points reduces the uncertainties in the 
fits, but does not break the basic redshift-temperature 
degeneracy for SED fits. 

The size of the uncertainties on the estimated param- 
eters L and Td for the Chapman et al. sources are high- 
lighted in Fig. 2 along with their redshifts: the uncer- 
tainties on the estimates for the other samples were dis- 
cussed by BBC03. The individual uncertainty on any 
determination is comparable to the extent of the 
parameter space spanned by the Chapman et al. sample. 
We also highlight the selection criteria for the SMGs in 
Fig. 2. The error bars to the Td are assigned by fitting 
SED templates to the colors of observed SMGs, using a 
maximum-likelihood estimator combined with broad pri- 
ors, and without imposing any condition on luminosity. 
The error bars in luminosity are derived from the range 
of luminosity values for which \ 2 < 4 for the two data 
points as the temperature moves away from the best- 
fitting value. Examples of correlations between fitted 
values of L and Td are shown by a joint simple x 2 fit to 
the flux density data for two typical SMGs in Fig. 3. 
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The different galaxy samples occupy somewhat differ- 
ent regions of the L-Td plane shown in Fig. 1 . Part of 
this separation is due to their distinct selection crite- 
ria. The low-redshift IRAS galaxies occupy a relatively 
well-defined strip, whether their SEDs are determined 
from either far-lR-submillimcter colors (open squares; 
Dunne et al. 2000) or from radio-far-lR template fit- 
ting (solid lines; Chapman et al. 2003a). The more dis- 
tant Stanford et al. (2000) sample includes galaxies with 
greater luminosities, and with typically higher and more 
scattered temperatures. These two samples match well 
at intermediate luminosities, L ~ 10 12 L Q , for which 
typical temperatures of about 40 K are determined. It 
is likely that both of these samples are biased against 
galaxies with cooler dust temperatures, because their in- 
clusion requires a detection at the relatively short ob- 
served wavelength of 60 /im. This corresponds to an even 
shorter restframe wavelength of about 45 /im for the typ- 
ical galaxy in the Stanford et al. sample. A few faint 
galaxies selected jointly at radio and mid-IR wavelengths 
are found to have SEDs consistent with lower tempera- 
tures (Garrett et al. 2001; see also Chapman et al. 2002). 
The new, large and more accurately defined Chapman et 
al. (2004a) SMG sample (roughly detected at a restframe 
wavelength of 250 /im) also is systematically cooler, with 
a locus of points that lie at significantly lower tempera- 
tures than the trend of the ZR^S-selected samples in the 
L-Td plane shown in Figsl & 2. 

The track of the L-Td distribution of SMG points in 
Fig. 2 appears to be relatively tight, but this does not 
necessarily reflect a tight intrinsic scatter in the proper- 
ties of the SMGs. Rather, the shape of this distribution 
is determined mainly by the condition of a detectable 
submillimcter-wave flux density, as shown by the dashed 
curves that trace the limits to detection for a 5-mJy 850- 
/im source in the Chapman et al. (2004a) SMG sample 
at a variety of redshifts in Fig. 2. Thus, it is impossi- 
ble to select SMGs in the high-Td-low-L region shown 
in the figures at any redshift in existing submillimeter- 
wave surveys. In a deeper survey, as will be possible with 
the Atacama Large Millimeter Array (ALMA; Wootten 
et al. 2001) the region of the parameter space that can 
be probed will expand, and many high-redshift galaxies 
will then probably be found in this region. The nature of 
high-redshift galaxies means that there are many possi- 
ble explanations for why they avoid the low-luminosity- 
high-tcmpcraturc region of the L-Td parameter space. 
They could have systematically different dust grain prop- 
erties and spatial distributions of dust as compared with 
their more luminous companions, or they could be dusty 
galaxies at the start or end of a burst of enhanced star- 
formation/ AGN activity, cither brightening from or fad- 
ing into the distribution of much more normal quiescent 
galaxies. Note that regular star-forming Lyman-break 
galaxies (LBGs; Steidel et al. 2004) will be found some- 
where in the L-Td plane at high redshifts (Fig. 2d). By 
detecting their lower far-IR luminosities using ALMA it 
should be possible to understand the currently ambigu- 
ous relationship between the LBGs and SMGs in more 
detail. 

The corresponding conditions for radio detection at 
1.4-GHz - that a galaxy must exceed a flux density of 
about 30/iJy is also shown in Fig. 2. The shape of the 
radio selection function is much less sensitive to tempera- 



ture, depending much more strongly on luminosity, with 
a more intuitive redshift dependence. 

The absence of points at high-luminosities in the L- 
Ti diagram reflects the steep high-luminosity cutoff of 
the luminosity function of dusty galaxies - there are few 
very luminous objects to find. Note that the observed 
strong luminosity evolution of dusty galaxies by a factor 
(1 + z)~< with 7 ~ 4 for z < 2 (Blain et al. 1999), is 
implicit in the results shown for different redshift ranges 
in Fig. 2. Very luminous galaxies are not represented in 
large numbers in the Stanford et al. sample at z ~ 0.5, 
yet are abundant in the SMG sample at z ~ 2.5. 

Chapman et al. (2004a) found wide discrepancies while 
comparing the scatter in the radio-submillimctcr flux ra- 
tios for SMGs with redshifts to the ratio expected as- 
suming the updated Carilli-Yun 1.4-GHz-850-/im red- 
shift estimator (Yun & Carilli 2002). The range of dust 
temperatures for SMGs shown in Fig. 2 is spread over 
a factor of at least 1.4 for SMGs with L ~ 1O 13 L . 
This implies a radio-submillimctcr far-IR photometric 
redshift accuracy no better than Az ~ 1 based on the 
temperature redshift degeneracy (BBC03) for the typi- 
cal SMG spectroscopic redshift z ~ 2.5: this is in marked 
contrast to recent claims of much more accurate photo- 
metric redshifts by Aretxaga et al. (2003). The pres- 
ence of hotter dusty galaxies at slightly lower luminosi- 
ties in the Stanford et al. (2000) sample (Fig. 1) hints 
that photometric redshift errors could be even larger 
if SMGs were selected down to deeper submillimeter 
flux density limits. Photometric redshift techniques at 
shorter near- and mid-IR wavelengths offer better hope 
of progress, either using redshifted polycyclic aromatic 
hydrocarbon (PAH)/silicate dust grain absorption and 
emission features, or stellar spectral features (Simpson 
& Eisenhardt 1999; Sawicki 2002; Efstathiou & Rowan- 
Robinson 2003). Regular multi-band near-IR, optical, 
UV photometric redshift techniques also hold promise, 
as many SMGs have the expected optical colors for their 
spectroscopic redshift (Chapman et al. 2004a), despite 
having far-IR luminosities that are almost impossible to 
predict from the optical data. Nevertheless, accurate po- 
sitions, perhaps from radio or ALMA observations, are 
still necessary in order to confirm that a particular faint 
galaxy with a photometric redshift really is the counter- 
part to a particular SMG. 

There is almost no overlap in physical properties be- 
tween well-studied low-redshift IRAS galaxies and the 
existing catalogs of SMGs, as shown in Figs 1 & 4. These 
classes of galaxies occupy totally disjoint regions of the 
L-Td parameter space, being separated by a factor of 
30% in temperature at comparable luminosities, and by 
a factor of 10 in luminosity at comparable temperatures. 
In part this is a selection effect, imposed by the require- 
ment of a 850-/im flux density in excess of 5 mJy at low 
redshifts and/or luminosities, as shown in Fig. 2. How- 
ever, the lack of overlap in the figure could also be a 
real astrophysical effect, perhaps reflecting either a more 
diffuse distribution of dust and stars or a lower typical in- 
terstellar medium metallicity in the high-redshift SMGs 
as compared with the IRAS galaxies. Spatially-resolved 
imaging of the radio, submillimeter or far-IR emission 
from these galaxies will be necessary to reveal the true 
relationship between these populations, its causes and 
the relative astrophysics of the IRAS and SMG galaxy 
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populations. 

We emphasize that spectroscopic redshifts for SMGs 
are essential to draw all these conclusions: without them 
reliable temperature and thus luminosity information 
cannot be derived. 

4. NEW POPULATIONS 

A significant fraction (~ 50%) of the SMG population 
with 850- /xm flux densities in excess of 5 mJy are included 
in the Chapman et al. (2004a) redshift sample (Fig. 2). 
Unless there is a significant systematic shift in the SEDs 
of SMGs at fainter flux densities (~ 1 mJy) , then slightly 
fainter, less luminous SMGs will typically populate the 
region of Fig. 1 that lies between the existing SMG sam- 
ples and the Dunne et al. (2000) low-redshift IRAS 
galaxies, providing a smooth distribution in luminosity 
between the extremely luminous high-redshift galaxies 
and more ordinary local dusty galaxies. These fainter 
SMGs at moderate and high redshifts will be hard to 
detect at submillimctcr wavelengths prior to the com- 
missioning of ALMA. 

High-redshift dusty galaxies with fainter 850-^m flux 
densities may have luminosities similar to those of known 
SMGs, but higher temperatures. These would be an im- 
portant class of objects, both very luminous and hotter 
than any population detected so far (Fig. 1). The IRAS 
survey is not deep enough to find such galaxies (Stan- 
ford et al. 2000), and small ISO samples do not appear 
to detect them in large numbers (Garrett et al. 2001; 
Franceschini et al. 2003). The deepest ISO 15- /im im- 
ages (Altieri et al. 1999; Metcalfe et al. 2003) may in- 
clude some examples, but redshifts, and supporting deep 
submillimcter and radio data are not common for this 
population. 

With T d ~ 60K and L ~ 10 13 L Q such hot galaxies 
at z ~ 2-3 could be detected in current deep radio sur- 
veys (Fig. 2), but not in existing submillimcter surveys. 
A factor of 10 increase in submillimcter sensitivity is re- 
quired to probe this region of the L-T^ diagram. Be- 
cause confusion noise sets the 850-/um sensitivity limit at 
about 2mJy in field galaxy surveys using current tele- 
scopes (Blain et al. 2002), higher angular resolution or a 
large gravitational lensing magnification, not just longer 
integrations, would be necessary. 

Optical spectra taken in parallel to Chapman et al.'s 
SMG redshift survey have provided redshifts for many 
optically faint galaxies that have weak radio detections 
but no submillimctcr emission (Chapman et al. 2004b). 
These galaxies' optical spectral properties are similar to 
those of high-redshift SMGs, without signs of power- 
ful AGN, so it is plausible to suggest that their radio 
emission is produced by star formation. These optically 
faint galaxies are excellent candidates for a hotter, high- 
redshift, dust-enshrouded galaxy population, which can- 
not be detected in submillimeter surveys owing to the 
galaxies' high dust temperatures. The surface densities 
of these galaxies are comparable to those of the SMGs, 
so a survey covering just a few tenths of a square degree 
should find a statistical sample of several tens. 

Such a population of galaxies would generate back- 
ground radiation intensity in addition to that from the 
known SMGs and their fainter counterparts. The extra 
intensity that they contribute would mostly increase the 
background intensity at wavelengths shorter than about 



200 fim, where the background reported by Schlegel, 
Finkbcincr & Davis (1998) is undcrpredicted both by 
models of the evolution of SMGs (Blain ct al. 1999, 2002), 
and by summing over the flux and surface densities re- 
ported for SMGs and their fainter dusty counterparts. 
Possible sources of 200- /im background radiation are dis- 
cussed by Blain & Phillips (2002): a 50% increase in the 
intensity of background radiation above that accounted 
for by extrapolating the known SMG population is pos- 
sible at about 100 ^m, and could be due to a hotter high- 
redshift population of dusty galaxies. 

The known selection effects for existing submillimctcr 
and far-IR surveys, observational hints at a spectroscopi- 
cally similar, but submillimctcr- faint, population of high- 
redshift radio source, and the intensity of the 100-200- 
/im extragalactic background radiation, are all consis- 
tent with the existence of a proposed population of hot 
high-redshift dusty galaxies. To be certain that we have 
revealed this population it is necessary to probe deeper 
into the luminosity function of dusty galaxies at high 
redshifts, an opportunity that is available most easily at 
shorter mid-IR wavelengths. 

5. SPITZER PROBES OF NEW POPULATIONS 

The recently launched Spitzer offers excellent sensitiv- 
ity and an unprecedentedly large field of view at mid 
and far-IR wavelengths. While it is still not absolutely 
clear how deep Spitzer can integrate before being limited 
by confusion noise, it will certainly produce much bet- 
ter images of the far-IR sky than any previous mission. 
Fig. 4 shows the region of the L-T^ plane that Spitzer 
will probe to the estimated confusion limits (Blain et al. 
2002) at various wavelengths. Deep Spitzer observations 
can thus investigate the whole of the unprobed region of 
the L-Td plane near 60 K and 10 13 Lq at 24-/im. In the 
unlikely event that a source confusion limit brighter than 
0.1 mJy prevents Spitzer from reaching this depth, then 
the larger 2.5-m aperture SOFIA 4 airborne observatory 
will certainly be able to probe this region of the L-T^ 
plane in very long exposures. 

The 24-/im flux density expected for a distant galaxy 
as a function of Td and L depends most strongly on the 
value of the a parameter that describes the mid-IR SED. 
However, increasing or decreasing a by 0.5 leads to less 
than a factor of two change in the predicted far-IR lu- 
minosity at 24- /xm expected from a galaxy at redshift 
z = 2.5 for Td > 40 K. If an extreme value of a = —2.9 
(representative of the steepest known mid-IR SED in 
Arp 220) is assumed, then only galaxy luminosities that 
are a factor of 5 times higher than the limits shown in 
Fig. 4 can be probed using a 24-/im observation to a cer- 
tain flux density limit. Hence, even for the most extreme 
SED, a Spitzer 24-/xm survey can still access the unex- 
plored region of the L-Td parameter space highlighted in 
Fig. 4 for galaxies at z ~ 2.5. 

At wavelengths longer than about 100 /im, confusion 
noise from faint undetected galaxies makes it difficult to 
probe the unexplored region at high redshifts unless the 
telescope used has an aperture greater than of order 3-m. 
Hence, to sample the full range of known high-redshift 
galaxy SEDs near their restframe emission peak will re- 
quire a next-generation 10-m-class far-IR telescope, such 

4 http://sofia.arc.nasa.gov 
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as the NASA SAFIR vision mission. 5 

Its sensitivity to less-luminous and hotter galaxies at 
lower redshifts will also ensure that Spitzer can be used 
to reveal whether the absence of low- luminosity, hot, rel- 
atively low-redshift galaxies in Fig. 4 reflects a real up- 
per limit to the luminosity-weighted dust temperature in 
galaxies, or whether it is a selection effect brought about 
by insufficient sensitivity in existing surveys. Note that 
because the line defining the 24-/im selection condition 
is predicted to be orthogonal to the track of the existing 
data points through Fig. 4, Spitzer observations will be 
sensitive to any high-redshift galaxies with luminosities 
in the rang e 10 n -10 12 L Q , and temperatures T d > 50K. 
Hence, they could readily test whether the upper tem- 
perature limit to the SEDs in existing samples is a real 
astrophysical result or a selection effect. 

A further benefit of the great sensitivity of Spitzer is 
its ability to detect any luminous galaxies that lie on 
the cold side of the SMG track in Fig. 4 (e.g. Chapman 
et al. 2002). The details of the detectability of cooler 
galaxies at the short rest wavelength of 24 /im are un- 
certain, and will remain so until the release of the first 
Spitzer data in mid-2004. Although such cold luminous 
galaxies should be detectable in existing submillimeter 
surveys, the much larger sky area that should be cov- 
ered by Spitzer - of the order of 100 deg 2 as compared 
with the 0.5-1 deg 2 covered by all existing deep mm and 
submillimeter-wave surveys - will yield a larger number 
of potential examples. Deep Spitzer images of SMG fields 
will also provide a direct test of the SEDs assumed for the 
SMGs based on the far-IR-radio correlation when plac- 
ing the SMGs on the L-T& plane in Fig. 1, and thus will 
test the evolution of this correlation. Note that to iden- 
tify the most interesting sources among the huge crop 
of Spitzer detections, deep radio and submillimeter data 
will still be very valuable, to provide accurate positions 
and SED/tcmperature limits respectively. 

If Spitzer detects faint high-redshift /iJy radio galaxies 
and confirms that they lie in the unexplored region of the 
L-T d plane (Fig. 4) , then we will have a full measure of 
the scatter in the SEDs/temperatures of high-redshift lu- 
minous dusty galaxies. If not, then we will have found a 
new population of distant radio-selected galaxies without 
obvious AGN signatures (Chapman et al. 2004b; Alexan- 

5 http://safir.gsfc.nasa.gov 



der et al. 2004). These galaxies would have to lie on the 
radio-loud side of the low-redshift radio-far-IR correla- 
tion. This new class of galaxies could be either X-ray and 
optically faint, dust-enshrouded, Compton-thick AGNs, 
or distant star-forming galaxies with anomalously pow- 
erful radio emission, which are found only at the 1% level 
in low-redshift samples (Yun et al. 2001). 

6. CONCLUSIONS 

The recent acquisition of a large sample of SMG red- 
shifts provides the opportunity to constrain the SEDs of 
the most luminous (> 10 12 Lq) high-redshift galaxies. It 
is now reasonably certain that the typical high-redshift 
SMGs selected at 850-/um are systematically cooler than 
local dusty galaxies of comparable luminosity, and are 
thus without direct low-redshift analogs. The scatter in 
their dust temperatures is about 50%, sufficient to in- 
troduce a ~ 30% — 40% uncertainty in any photomet- 
ric redshifts derived from far-IR, submillimeter, or radio 
photometry, no matter how accurate. The true scatter 
of the dust temperatures in high-redshift dusty galaxies 
certainly could be still greater, but hotter examples can- 
not be selected by existing submillimeter surveys. Faint 
radio galaxies that are undetected in the submillimeter, 
but whose optical properties are like those of the SMGs 
hint at the existence of a hotter population, with lumi- 
nosities L ~ 10 13 Lq and temperatures T<j — 60 K. These 
galaxies would produce a luminosity density compara- 
ble to that of the existing populations of SMGs. The 
deepest forthcoming mid-IR Spitzer surveys should soon 
confirm whether this population exists, while providing 
a complete inventory of the range of SEDs present in 
high-redshift dusty galaxies in fields that have excellent 
supporting long-wavelength data. 
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Fig. 1. — Range of luminosity L and temperature values determined for well-studied dusty galaxies at a variety of redshifts. Open 
squares, low-redshift IRAS galaxies (Dunne ct al. 2000); open circles, modcrate-rcdshift very luminous IRAS galaxies (Stanford ct al. 
2000); open triangles, high-redshift ISO galaxies (Garrett et al. 2001); filled circles, spectroscopically identified SMGs (Chapman et al. 
2003b, 2004a). Crosses and plus signs show galaxies with well-determined SEDs from a wide variety of sources (including SCUBA and 
ISO surveys). Crosses represent a selection of gravitational lenses, while plus signs represent unlcnscd galaxies (for both see BBC03); and 
spectroscopically identified SMGs - filled circles (Chapman ct al. 2003, 2004a). The sample represented by diagonal crosses includes M82, 
the Milky Way, NGC958 and Arp220, which are shown by the largest crosses between 2 X 10 10 Lq and 1.6 X 10 12 Lq respectively. For the 
non-Chapman samples larger symbols reflect more accurate determinations: sec BBC03 for a description of the errors for these samples. See 
Fig. 2 for uncertainties in the estimates for the Chapman et al. sample, for which all symbols are the same size. The thick solid line shows 
the result derived by Chapman et al. (2003a) for the L-T<j values of low-redshift IRAS galaxies. The interquartile range for this sample is 
bracketed by the two thinner solid lines, the separation of which agrees well with the range of properties of the Dunne et al. sample. The 
dashed and dotted lines show the fits obtained to merging and quiescent low-redshift IRAS galaxy data respectively by Barnard (2002). 
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Fig. 2. — SED parameters for the galaxies shown in Fig. 1 divided by redshift range, and the errors on the parameters fitted to the 
Chapman ct al. (2003b, 2004a) sample, symbols arc as in Fig. 1. The selection effects that limit the SED parameters for galaxies in this 
sample arc represented by the lines at submillimctcr and radio wavelengths for the labeled rcdshifts. Galaxies can be detected to the 
high-luminosity side of each line. There is a trend for the most luminous galaxies to be at the highest redshifts, reflecting the strong 
evolution of the luminosity function of the dust-enshrouded galaxy population. Nevertheless, note that broad ranges of luminosity, and 
more importantly temperature, are sampled in each redshift interval regardless of the selection effects. 
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Fig. 3. — Examples of the uncertainties in the SED fit obtained for two typical radio-submillimctcr galaxies in the HDF field. Left: 
SMMJ123707.21+621408.1, No. 242; right: right is SMMJf2372f. 8+621035.3, No. 378) from Chapman et al. (2004a), at z = 2.48 and 
1.51 respectively). The error bars match the results plotted in Fig. 2, and are determined by making a full maximum-likelihood fit to the 
temperature for an SED of arbitrary luminosity and then finding the range of best-fitting values of L that have a value of \ 2 < 4 at that 
value of Td. The contours show reduced x 2 values of 1, 4, 9, 16 and 25 derived directly for comparison. The values are consistent, 

but slightly distinct, reflecting the different fitting procedures. The minimum x 2 value is denoted by the star. Note that the extent of the 
contours in the temperature direction can be reduced by adding any data from Spitzer mid- or far-IR observations. 
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Fig. 4. — Full range of L-T^ values determined for dusty galaxies (see Fig. 1). The lines show the limits for selecting galaxies in surveys 
close to the confusion limit for a 0.85-m aperture telescope like Spitzer and a 3-m-class telescope like Herschel or SOFIA, for a galaxy at 
rcdshift z = 2.5, typical of the identified SMGs. The solid lines show the limits of Spitzer at 24 fim (thick lines; 0.1 mjy), 70 fim (medium 
lines; 5mjy) and 160 ^m (thin lines; 80mjy). The dashed lines represent the limits of a 3-m aperture telescope at 250 fim (thin lines; 
70mjy), and 100 fim (thick lines; 1.5 mjy). The selection limits in the existing Chapman et al. samples (Fig. 2) are shown by the dotted 
lines. Galaxies to the high-L side of all the lines can be detected. Reaching the unexplored region to find ultraluminous (L ~ 10 Lq) 
galaxies at ~ 60K requires Spitzer observations at 24 fim or an aperture larger than about 3 m at observing wavelengths longer than 
100 fim. To probe a significant part of the unexplored region of the diagram at submillimeter wavelengths would require a survey depth of 
order 0.1 mjy at 850 fim with ALMA. 



